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ABSTRACT 

^ ■ In recent years, the Extremely Low Mass White Dwarf (ELM WD) survey has 

p f quintupled the number of known close, detached double WD binaries (DWD). 

The tightest such DWD, SDSS J065133.33+284423.3 (J0651), harbors a He WD 
eclipsing a C/0 WD every ~ 12min. The orbital decay of this source was 
recently measured to be consistent with general relativistic (GR) radiation. Here 
we investigate the role of dynamic tides in a J0651-Like binary and we uncover the 
potentially new phenomenon of "antiresonance" locking. In the most probable 
scenario of an asynchronous binary at birth, we find that dynamic tides play a 
key role in explaining the measured GR-driven orbital decay, as they lock the 
system at stable antiresonances with the star's eigenfrequencies. We show how 
such locking is naturally achieved and how, while locked at an antiresonance, 
GR drives the evolution of the orbital separation, while dynamic tides act to 
synchronize the spin of the He WD with the companion's orbital motion, but only 
on the GR timescale. Given the relevant orbital and spin evolution timescales, 
the system is clearly on its way to synchronism, if not already synchronized. 

Subject headings: binaries: eclipsing — binaries: general — stars: interiors — 
stars: oscillations — white dwarfs — gravitational waves 



1. Introduction 



Short period DWDs are the most numerous sources for the next ge neration of low- 

frequency (mHz) sp ace-based gravitat ional wave (GW) detectors (e.g. LISA. lDanzmann fc the LISA study 
19961 . iHughe jbood . and eLISA/NGO. lAmaro-Seoane et al.l[2012l l. The eclipsing DWD SDSS 
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J065133.33+284423.3 fhe reafter J0651. ISrown et al.lboill: iHermes et alJ l2012h. discovered 



2012 



Kilic et al. 



2010 



201 ll ). harbors a tidally 



by the ELM WD survey flBrown et aLll2010l 
deformed He WD orbiting a C/0 WD every ~ 12min; this is the shortest-period detached 
system known. Here, we investigate the effects of dynamic tides on the orbital evolution of 
this binary. 

Previous investigations of tides in J0651 placed limit on the compon ents' tidal p arameter 
Q, assuming tidal heating to be the cause of the observed luminosity ( lPirdl201ll ). or used 
polytropes to estimate the deviation fro m the pure GR-dri ven inspiral, under the assumption 
that the components are synchronized (IBenacquistal 1201 ll ) . 



T he orbital evolu t ion of detached DWDs is often assumed to be driven by GR emission 
alone. IWillems et al.l (120 lOf ) validated this assumption provided that the binary is in the 
regime of dissipative quasi-static tides (in the limit of orbital and rotational periods ^ 
y/W{GM)^). The inefficiency of quasi-static tides means GR emission eventually drives the 
system out of synchronism and DWDs enter a regime in which dynamic tides and resonantly 
excited gravity modes ((^-modes) must be taken into account. 

Here, we use our newly developed tidally excited stellar pulsation code CAFein (Valsec- 
chi et al. to be submitted) and detailed models of a J0651-Like He WD to study in detail 
the effect of radiative dissipation of dynamic tides on the orbital evolution of J0651. Our 
focus on radiative dissipation follows our u nderstanding of dissipative dynamic tides in non- 
degenerate stars with radiative envelopes (IZahnlll975l : IWitte &: Savonijd 120011 ). 



Fuller fc Lail (120111 |2012| ) investigated the impac t of resonant l y exc ited (7— modes on 



the orbital evolution of C/0 DWDs. Dissipation in Fuller fc Lail ( 12012! ) was treated via 
the so-called outgoing wave boundary condition (BC), which implicitly assumes that tides 
are efficiently damped via radiative damping or non-linear effects. The authors focused 
on the evolution of (7— mode amplitudes during resonance, without a detailed treatment of 
tidal dissipation in and out of resonance, and found that resonant mode excitation could 
quickly drive the binary to (near) synchronization of the spin and orbit on the dynamic- 
tides timescale. Here, we reach fundamentally different conclusions: we find that, for a 
broad range of initial spin frequencies, dynamic tides act to drive the spin-orbit frequency 
mismatch to and trap it in an aniiresonance, where the dissipation due to tides effectively 
vanishes. The orbit subsequently evolves due to GR emission, and the only action of tides 
is to keep the spin-orbit frequency difference trapped at the antiresonance value; as GR 
emission proceeds, both the spin and orbital frequencies increase, while their difference is 
locked, and the system is driven to synchronism, but this process takes place on the GR 
timescale. 
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Very recently, iHermes et al.l ( 120121 ) reported the detection of orbital decay in J0651 
consistent with the GR prediction. Here we show that dynamic tides potentially play a key 
role in such measurement, as they naturally drive the system into antiresonance locking, a 
new phenomenon we identify, with further evolution of the orbit occurring only on the GR 
timescale. 



2. Observed and Model's Properties 



The DWD system J0651 was originally discovered by lBrown et al.l ( 120111 ). iHermes et al. 



( 120121 ) measured an effective temperature, gravity, radius, and mass for the He WD of 
Teff = 16530 ± 200 K, \og{g) = 6.76 ± 0.04 dex, Ri = 0.0371 ± 0.0012 i?0, and Mi = 
0.26 ± 0.04 Mq; they also measured the C/0 WD mass to be 0.5 ± 0.04 Mq and the orbit 
is circular. The measured orbital period and orbital decay are P = 765.206543(55) s and 



P 



-9.8 ± 2.1 



X 



10 ss , respectively. For comparison, the orbital period decay 
ex pected from GR emi ssion alone is P = (—8.2 ± 1.7) x 1 0~^^ss~^ for the masses derived 



bv IHermes et al.l (l2012l ). We use MESA ( IPaxton et al.ll201ll ) to create a He WD model with 
Teff ~ 16766 K, log(^) = 6.68 dex, i?i = 0.0363 i?©, and Mi = 0.23 M©. These parameters 
are within 1 a of the measured Mi and Ri and 2 cr of the measured Teff and \og{g). 



3. Tidal Evolution 

3.1. Non-Adiabatic Stellar and Tidal Pulsations 

When tides are dynamic, the tidal forcing exerted by a star in a binary can interact 
with one or more of the companion's eigenfrequencies. Here we take the He WD (primary) 
to rotate uniformly with an angular velocity Qi and, given the MESA model's current cool- 
ing timescale tcooi (see below), we assume that its stellar parameters do not evolve. We 
also neglect the effects of Coriolis and centrifugal forces (the break-up spin for the model 
is ~ 2.5 min). Under these assumptions, the star's eigenmodes correspond to those of a 
non-rotating spherically symmetric star in hydrostatic equilibrium. We treat the C/0 WD 
(secondary) as a point mass. Let r = (r, 6, 0) be a system of spherical coordinates with 
respect to an orthogonal frame co-rotating with the primary. The system of eq uations gov- 



erning non-adiabatic and non-radial stellar pulsations has been derived e.g. by lUnno et al. 



( 1l989l ) (Eqs. 24.7 - 24.12 in §24). It comprises six homogeneous ordinary differential equa- 



tions involving the complex dimensionless eigenfrequency oo {a = u^jGM/ R^, where a is 
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the dimensional eigenfrequency) . It is convenient to use 

yi = —, y2 = — [- + ^ ] = , z/3 = 

r 9^ \P J 9 ^ 9f 

yi = 1/5 = — , 1/6 = -j— (1) 

9 dr Cp Lr 

as integration variables in the eigenvalue problem. Here C,r and are the radial and orthogo- 
nal component of the displacement of the star's mass element from the equilibrium position, 
p the density, Cp the specific heat at constant pressure, and p', 6S, and 5Lr the perturbed 
pressure, potential of self-gravity, entropy and radiative luminosity, respectively. The effect 
of convection on the oscillations is neglected. 

The central BCs are that the entropy is constant during a single oscillation {6S = 0) 
and that {p'/p + $'), and are regular. The surface BCs are determined by assuming 
that near the surface the pressure and the density drop steeply outward, by considering that 
there is no inward radiative flux, and by requiring the continuity of $' and its first derivative. 

The tidal action from the companion is included as a small time-dependent force which 
perturbs the hydrostatic equilibrium of a spherically symmetric static star. The system 
of equations describing the tidal response of the star is formally identical to the equations 
describing pure stellar pulsations, provided that the perturbed potential of self-gravity $' 
is replaced by the total perturbation of the potential = $' + ex W^, where erW is the 



tide-g enerating potential. We expand in a Fourier series as (see e.g. iPolfiiet &: Smeyers 



19901) 



eTW{r, t) = -ex ^ ^ ^ Q,„,fc f — j Yl^^O, 0)exp[i(cT„,fct - kQorhr)] 

1=2 m=-lk—~- ^ 



(2) 



where Yi"^{6, 0) are unnormalized spherical harmonics of harmonic degree /, and azimuthal 
number m, and k is the Fourier index; 6 is the colatitude in the WD's co-rotating frame, and 
is the longitude. At time t = 0, the angle = marks the position of the periastron. The 
dimensionless parameter ex = (-Ri/a)^(M2/Mi) gives the ratio of the tidal force to gravity 
at the star's equator. The semi-major axis of the binary is a, QqA is the mean motion, r 
the time at periastron passage, am,k = ^^orb + "^^i is the tidal forcing angular frequency 
with respect to the co-rotating frame, and Q^m,fc are Fourier coefficients. For a binary with a 
circular orbit (like J0651) the only non-zero Q^m,fc are the ones with k = — m and are defined 
as 

1-2 



c,„=|^pr(o)(^)" (3) 

(/ + \m\y. V « / 
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where Pl^{x) are Legendre polynomials of the first kind. Refer to lPolfliet fc SmeyersI (119901 ) 
for further discussion of Ci^m,k- The expansion ([2]) shows that the tidal action from the 
secondary induces in the primary an infinite number of forcing angular frequencies (Jm^k- 
Tides are dynamic when am,k 7^ 0. 

In principle, all values of /, m, and k must be considered in the set of equations describing 
tidally excited stellar pulsations. However, as Q^m,fc (-Ri/^)'"^) investigations on dynamic 
tides are often restricted to the terms belonging to / = 2. Here we consider only the leading 
term in the expansion of e^W: {l,m, k) = (2,-2,2). 



3.2. Secular Evolution of the Orbital Separation and Stellar Spin 



Here we summarize the main equations describing the secular evolution of a and Qi 
du e to dynamic t ides for a b inary with a circular orbit and for {l,m,k) = (2,-2,2); refer 



to IWillems et al.l (120031 . |2010| ) for a more detailed derivation and for an explanation of the 
various parameters entering such equations. 

The primary's tidal deformation perturbs the external gravitational field, and thus the 
Keplerian motion of the binary components . In the framework of the theory of osculatiri g 
elements in celestial mechanics (jSternelll960l : iBrouwer fc Clemencelll96ll : iFitzpatrickl Il970h . 
the rate of change of the orbital semi-major axis due to the tidal deformation of the star is 
given by 



/ da 
\di 



Svr M2 



where G^l(0) 



-2mP^I(0)Q, 



a 

for e 



/+3 



\Fi^rn,k\^iT^'yi,m,kG\'^l^,^{0) (4) 

0, Porb is the orbital period, and the dimen- 



sionless -F),m,fc are defined as 

PLm,k 



Rl '^l,m,k{Rl 



+ 1 



I jn.k I ^ 



(5) 



GMi eTCi^rn,k 

and measure the response of the star to the various tidal forcing frequencies; Fi^m,k can be 
computed from a solution to the modified pulsation equations described above. The angle 
'li,m,k describes the phase-lag between the tidal forcing and the response of the total potential 
induced by the dissipation in the system. 

This phase shift produces a torque on the tidally deformed star, which acts on its spin 
yielding 



'd^i 



Stt / a Ml Ml \ Ma a^l^ ( Ri 



Porb + 



Ml h 



/+3 



— P/ 



(6) 
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4. The Dynamic tides regime in a J0651-Like binary 



10.000 



1.000 r 



0.100 



0.010 



0.001 




Fig. 1. — Left y-axis: radial profile of the Briint-Vaisala (A^) and Lamb (L2) frequen- 
cies squared (black-solid lines). The dots mark the positions of the nodes for the modes 
9i, 95, 9w, 9i5, 920, 930, 9m, 960, 9so- The horizontal grey lines mark the position of the tidal 
forcing frequencies for spins (from top to bottom) Qi = 0.15 QqA, ^1 = 0.5 QoA, and Qi 
= 0.9 fiorb- All the frequencies are normalized by a/ GM/R^. Right y-axis: the dashed-dot 
line shows the ratio of the star's thermal to dynamical (rdyn/Tth) timescales. Where Tdyn/fth 
is small non-adiabaticity becomes relevant. 



Here we first examine whether dynamic tides are relevant in a J0651-Like binary by 
comparing the He WD eigenfrequencies with the tidal forcing frequencies. We use CAFein 
(Valsecchi et al. to be submitted) to compute the non-adiabatic / = 2 eigenfrequencies for 
(7 — modes between 9i-9%o. The complex eigenfrequencies are summarized in Tabled! In 
Fig. [1] we show the propagation diagram for a J0651 -Like He WD, where the squared Briint- 
Vaisala frequency is A^^ = g{V^^ dhip / dr — d\n.p / dr) (IBrassard et al.lll99ll ). while the squared 
Lamb frequency is given by = l{l + l)Cs/r^. The (7-modes propagate in regions where 
Lo'^ < N"^ and Lf. The dots in Fig. [T]mark the positions of the radial nodes for the modes 
considered in Table (H while the horizontal grey lines mark the position of the tidal forcing 
frequency (w^) if ^1 is varied. As and cover the same frequency regime, resonances 
between tides and eigenmodes can be relevant for a J0651-Like binary, and hence tides are 
dynamic in this system. 
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Table 1: Non-adiabatic eigenfrequencies for a J0651-Like He WD. The subscripts "R" and 
"I" denote the real and imaginary part, respectively. Note that |u;i| increases with the order 
of the mode, indicating that non-adiabaticity becomes relevant for high-order modes (see 
Fig. d]). Also, ui < for modes up to gig, indicating that these modes are unstable. 



mode 








8.928 


10- 


-1 


~ 10^14 




95 


3.915 


10- 


-1 


-1.65-10- 


ll 


9io 


2.379 


10- 


-1 


-1.057-10 


-9 


915 


1.701 


10- 


-1 


-7.804-10 


~9 


920 


1.326 


10- 


-1 


8.403-10- 


-9 


930 


9.237 


10- 


-2 


1.959-10- 


'6 


940 


7.065 


10- 


-2 


9.882-10- 


-6 


960 


4.785 


10- 


-2 


2.842-10- 


-5 


980 


3.611 


10- 


-2 


5.519- lo- 


-5 



5. Synchronism Through Dynamic Tides : the Antiresonance Locking 

In Fig. [2] we show the timescales for the evolution of a and Qi due to tides calculated 
using Eqs. (jlj) and We fix QoA and M2 to the currently observed values and vary Qi, thus 
changing ut- Here we consider Qi ranging between ~ (0.15-0.999) fiorb- Black and grey data 
points represent positive and negative timescales, respectively. As expected, dgec < when 
^i,sec > and vice versa. Note that Qi evolves faster than a by about 2 orders of magnitude 
. Dips in the timescales correspond to resonances between the tidal forcing frequencies and 
the star's eigenfrequencies, while the peaks mark the "antiresonances" . 

The sign of the timescales depends on the angle that parametrizes tidal dissipation, 
li,m,k- This angle, in turn, is related to the physical angle describing the orientation of the 
bulge in the star's reference frame via 

</> = -li,m,k/m, (7) 

or for m =-2, = 7/2. 

When 72,-2,2 = or TT, the rate of change of a and Qi becomes zero and the timescales 
display an antiresonance. At these lag angles the bulge is either aligned with the companion 
or 7r/2 out of alignment; in either case there is no tidal torque and therefore no change in the 
orbit or spin due to tides. On the other hand, resonances occur when 72,-2,2 equals either 
7r/2 or — 7r/2, and thus (p equal to 7r/4 or — 7r/4; in these configurations, the torque is at a 
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Fig. 2. — Timescales for the secular evolution of orbital separation and stellar spin due 
to dynamic tides. Top: log\ta,tid\ = log|a/asec|- The horizontal dashed line shows the 



theoretically expected GR timescale log|ta,GR| from iPeterd (119641 ) for the masses considered 
here (ta,GR — 5Myr). Bottom: loglt^J = log|f2i/f2i sed- Grey and black lines and data 
points correspond to negative and positive timescales, respectively. As a reference, the first 
dip on the right marks the gi—mode resonance. Note that the x-axis is not uniform, and 
that I to. 



10-2 \t 



a,tid • 



maximum. 



The agreement between the observed orbital decay and the GR prediction implies that 
tides are not currently affecting the orbit. If the He WD were born close to synchronization or 
evolved to synchronism with the companion's orbital motion, tides are quasi-static (cut ^ 1) 
or static (cut = 0) and the orbital evolution tim escales due to quasi- static tides are longer 



than the theoretical GR timescale (see Fig. [2] or IWillems et al.ll2010[ ). On the other hand, 
if the He WD were born in or evolved to an asynchronous state, the signs of the timescales 
in Fig. m can result in spin evolution that stably locks the star at antiresonances, where 
^a,tid > ^a,GR (see details below). The crossing of the antiresonance peaks with the horizontal 
dashed line in Fig. [2] is visible only for cut — 0.23 because of the resolution adopted during 
the scan of the parameter space in Qi, but ta.tid ^ i^a.cR whenever the angle ■yi^m,k in Eq. (jl]) 
passes though and vr. 

Fig. E] describes in detail the spin evolution that results in antiresonance locking. For a 
binary on the left of the leftmost peak in the figure, where Qi < 0, dynamic tides act on Qi 
and move the binary to the right eventually reaching the antiresonance. At antiresonance 
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Fig. 3. — Resonances with the modes gg and qiq and antiresonance locking. The various 
symbols are as in Fig. [2] with the addition of 72,-2,2 and a cartoon depicting the orientation 
of the tidal bulges in the star's frame (the black- filled circle denotes the companion). The 
arrows indicate the directions towards which dynamic tides drive a J0651-Like binary (see 
text); antiresonances where < and Hi > on the left side and right side, respectively, 
are trapping. Once trapped, ut = 2 (fiorb — ^1) is constant, while the subsequent GW 
emission implies fiorb = f^i > 0, so that Qi/Qorh 1- Note that the x-axis in this figure is 
not uniform. 



^a,tid > ^a,GR and GR drives orbital shrinkage without affecting Qi. As Qorh increases and 
the binary moves to the right of the peak, dynamic tides act on the spin to bring it back at 
antiresonance. This stable antiresonance locking occurs for each antiresonance with f2i < 
and f^i > on the left side and right side, respectively. On the other hand, the antiresonance 
in the middle panel of Fig. [3] and all regular resonances are unstable, as dynamic tides quickly 
drive the binary out of them. 



Mathemat ically, antiresonance loc king is analogous to typical resonance locking (as dis- 
cussed e.g., by IWitte fc Savonijd 119991 ). If the star is locked at an antiresonance frequency 



- 10 - 



Ua, then = Ua- Assuming the background model doesn't evolve yields = and thus 
CUT = 0. The latter assumption is justified here as for the MESA model tcooi — 90 Myr and 
it is longer than the relevant timescales considered. As cjt is constant, Cl^rh = and, since 
at antiresonance locking GR shrinks the orbit, fiorb = f^i > 0. The latter implies that the 
combination of orbital shrinkage due to GR and of antiresonance locking due to dynamic 
tides brings the He WD to synchronization with the companion's orbital motion on the GR 
timescale. In fact, defining r = f2i/f2orb, "we find f oc fiorb(^orb — ^i) > 0. Since the absolute 
difference between Qi and fiorb is fixed by the antiresonance condition and both and 
|r2orb| are increasing, r — )■ 1 on a timescale tj- = r/r = — (2/3)tGR^i(^orb — ^i) 

Consequently, the binary, while locked in antiresonance, it is driven to synchronism 
because of the GR-driven orbital decay. Whether J0561 has achieved synchronism already 
or it remains locked at antiresonance, still on its way to synchronism depends on the relative 
magnitude of the timescales involved (tr, tcR, and tcooi), as well as their history, linked to the 
initial properties of the system. Answering this question requires coupled time intergration 
for a range of initial conditions, an analysis that is beyond the scope of the present study 
focusing on the new phenomenon of antiresonance locking; such a study will be the subject 
of future work. 



u s s u 




Fig. 4. — Possible combinations for 72,-2,2 in crossing two consecutive resonances (R), when 
the crossing occurs at 72,-2,2 = — 7r/2. Black and grey lines denote positive and negative 
72,-2,2, respectively. We use "S" and "U" to mark stable and unstable antiresonances. 

Even though it is difficult to identify all the stable antiresonances in Fig. [2] at the 
shown resolution, the number of antiresonances expected between two resonances can be 
determined from the change in the angle 72,-2,2 in the passage through two consecutive 
resonances, as shown in Fig. HI This illustration shows all possible combinations in 72,-2,2 
when the resonance occurs at 72,-2,2 = — vr/2 (a similar picture can be drawn for resonances 
occurring at 72,-2,2 = 7r/2). The number of antiresonances expected is either none, two, 
or four. For the latter two cases, we expect alternating stable and unstable antiresonances. 
The effect of resonances and unstable antiresonances is to drive the binary towards the 



closest stable antiresonace. This behavior is generic for any object under the influence of 
non-adiabadic dynamic tides, if the evolution of its properties are neglected; the locations of 
the resonances and antiresonances depend on the details of the object's structure, but the 
general picture of antiresonance locking depends only on the continuity of the phase 7;,™,,^ 
between resonances. 



6. Conclusions 

In this Letter, we show that dissipative dynamic ti des potentially play a key role in 



explaining the recently measured orbital decay of J0651 (iHermes et al.ll2012l ). even though 
their effects are not directly apparent. In fact, in the most probable scenario in which 
the binary was asynchronous at birth, the agreement between the observed orbital decay 
and the GR prediction is achieved through the new phenomenon of antiresonance locking 
identifled here. Assuming the He WD doesn't evolve on a relevant timescale, dynamic tides 
naturally lock the system at stable antiresonances with the star's eigenfrequencies. While 
locked, GR drives the evolution of the orbital separation, with dynamic tides maintaining 
the antiresonance condition and drive the binary to synchronism on the GR timescale. Time 
integration of the binary orbit for a wide range of initial conditions and a self-consistent 
calculation of all the relevant timescales will allow to assess whether the system is still locked 
and away from synchronism or if it has had time to reach synchronism. Future measurements 
of the He WD spin can potentially provide a check of the theoretical predictions presented 
here. 
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